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AbSMCt 

The  dielectric  properties  of  and  dry  natural  snow  were  studied  In  the 
trequency  range  of  50  Hz  to  100  kHz  to  determine  whether  measurements 
made  in  this  frequency  range  might  prove  useful  in  evaluating  the  water 
content  of  snow.  Dieiectric  heating  at  20  kHz  proved  a  very  useful  means  of 
modifying  the  water  content  from  0  to  30%  by  w^ght.  Six  dittorent  natural 
snows  were  used  in  these  experiments.  Mettwotn^  was  onoiyzed  tor  conduc¬ 
tivity,  pH,  and  impurity  content.  In  addition  to  developing  information  on  the 
dielectric  properties  of  w^  and  dry  snow,  we  measured  the  changes  produced 
in  dry  snow  by  altering  Its  density  overthe  range  ofO.1 1  to0.66g/cm^.  D^ils 
of  the  experimental  technique  and  the  data  obtained  are  tolly  reported.  Our 
results  do  not  lead  to  optimism  about  the  usetolness  of  measurements  in  this 
frequency  range  alone  tor  the  determination  of  water  content. 


For  conversion  of  Si  metric  units  to  U.STBritlsh  customary  units  of  measure¬ 
ment  consult  ASTM  Standard  E380,  Metric  PmcHce  Guide,  published  by  the 
ArTrertoon  Society  for  Testing  and  Materials,  1916RacsSf.,  Philadelphia,  Po. 
19103. 
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Determination  of  the  Water  Content  of  Snow  by 
Dielectric  Measurements 

PAUL  R.  CAMP  AND  DAVID  R.  LABRECQUE 


INTRODUCTION 

Knowledge  of  the  free  water  content  of  snowfields  is 
important  in  snow  mechanics,  hydrology,  avalanche 
control  and  in  satellite  imagery  interpretation.  Various 
techniques  for  determining  the  water  content  of  snow 
have  been  employed,  but  these  measurements  must  be 
made  by  teams  in  the  field.  Direct  methods  reviewed  by 
Boyne  and  Fisk  (1990)  include  melting  and  freezing 
calorimetry,  alcohol  calorimetry  and  ionic  dilution. 
Electrical  methods  have  focused  on  measurements  of 
the  dielectric  properties  at  frequencies  at  which  the  di¬ 
electric  constant  of  water  is  high  and  that  of  ice  is  low, 
generally  above  1  MHz.  These  methods,  reviewed  by 
Denoth  et  al.  (1984),  also  involve  field  studies  because 
they  require  an  independent  measurement  of  density. 
Qearly,  it  would  be  highly  desirable  to  be  able  to  make 
free  water  determinations  by  remote  sensing  or  by  auto¬ 
matic  instruments  that  can  be  interrogated  remotely 
without  the  need  for  field  measurements. 

A  wide  spectrum  of  frequencies  has  been  explored 
for  electrical  methods  of  determining  thfc  free  water 
content  of  snow,  the  most  popular  being  the  UHF  or 
microwave  ranges.  It  is  here  that  water  has  a  high  dielec¬ 
tric  constant  (about  88  at  0°C)  and  ice  a  low  dielectric 
constant  (about  4),  thus  providing  good  contrast  to 
reveal  the  presence  of  water.  Typically  these  measure¬ 
ments  focus  on  the  behavior  at  one  frequency  only.  In 
this  study,  we  have  chosen  to  see  what  might  be  learned 
from  measuring  the  dielectric  response  over  the  low- 
frequency  range  in  which  ice  has  a  Debye  dispersion. 
Our  thought  was  that  by  knowing  the  dispersion  prop¬ 
erties  of  ice,  measurements  at  several  frequencies  over 
the  dispersion  region  might  provide  simultaneous  equa¬ 
tions  the  solution  of  which  would  yield  the  water  con¬ 
tent  without  an  independent  measurement  of  density. 
The  argument  is  given  in  more  detail  below. 


Figure  1  shows  the  approximate  dielectric  constant 
for  ice  at  0°C  as  a  function  of  frequency  (solid  line).  The 
dielectric  constant  of  water  is  also  shown  (dashed  line). 
It  would  appear  that  one  might  easily  determine  the 
water  content  of  the  snow  by  combining  a  very  low  fre¬ 
quency  measurement  of  effective  dielectric  constant  (in 
which  the  contributions  of  the  water  and  the  ice  are  near¬ 
ly  equal)  and  a  high-frequency  measurement  in  which 
they  are  very  different.  However,  a  variety  of  things  dis¬ 
tort  the  low-frequency  end  of  the  ice  spectrum,  as  is 
shown  by  the  doned  line.  These  include  impurities  and 
electrode  polarization  effects.  As  frequency  increases, 
the  importance  of  these  anomalies  decreases  but  we  get 
into  the  dispersion  part  of  the  curve  where  the  calcula¬ 
tions  b^ome  more  difficult. 

Figure  2  shows  the  approximate  dispersion  with  fre¬ 
quency  of  the  conductivity  of  pure  ice  (solid  line)  and  of 
pure  water  calculated  from  the  relaxation  (dashed  line). 
Again  one  might  hope  to  be  able  to  infer  the  relative 
amounts  of  water  and  ice  in  a  sample  by  measurements 
at  two  fiequencies,  one  well  below  the  ice  dispersion 
and  one  at)Ove.  However,  as  is  shown  in  Figure  2,  the 
low-frequency  end  of  the  dispersion  flattens  out  due  to 
the  DC  conductivity  of  ice  and  to  the  presence  of  other 
low-frequency  dispersions.  Moreover,  both  the  DC 
conductivity  of  ice  and  the  strength  of  the  dispersion  are 
sensitive  to  impurities  in  the  ice.  More  important  still, 
the  DC  conductivity  of  water  is  dominated  by  ions.  The 
thermal  dissociation  of  water  alone  gives  pure  water  a 
resistivity  at  room  temperature  of  5  x  10^  ohm-m.  Dis¬ 
solved  CO2  fi’om  the  atmosphere  can  decrease  this  to 
2  X  10^  ohm-m.  The  respective  conductivities  are  also 
shown  in  Figure  2  as  dashed  lines.  Since  additional  im¬ 
purity  will  raise  the  conductivity  even  further,  this 
baseline  is  very  strongly  dependent  on  impurities.  We 
conclude  that  the  measurement  of  dielectric  constant  is 
more  promising  than  that  of  dielectric  loss. 


Frequency  (Hz) 

Figure  1.  Approximate  dielectric  constants  of  ice  and  water  as  a  Junction  of 
frequency  at  OX!  (extrapolated from  Camp  et  al.  1967). 


Frequency  (Hz) 

Figure  2.  Approximate  electrical  conductivities  of  ice  and  water  as  a  function  of 
frequency.  Dashed  line  to  the  right  shows  the  conductivity  of  pure  water  cal¬ 
culated  from  its  relaxation  spectrum  ( use  scale  at  the  right  for  these  values )  (data 
from  Camp  et  al.  1967). 
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Although  the  low-fre¬ 
quency  measurements 
may  be  distorted  by  elec¬ 
trode  and  other  phenom¬ 
ena  as  mentioned,  it  now 
appears  possible  to  use 
measurements  in  the  dis¬ 
persion  region  in  a  fairly 
simple  way.  The  advent 
of  inexpensive,  hand¬ 
held,  programmable  com¬ 
puters  makes  possible  the 
rapid  interpretation  of  this 
more  complicated  data. 
It  can  be  automated  in  the 


field,  almost  in  real  time. 

As  shown  in  Figures  1  and  2,  the  dielectric  properties 
of  ice  change  rapidly  in  the  dispersion  range,  while 
those  of  water  are  essentially  constant.  Thus  we  should 
be  able  to  separate  the  effects  of  water  and  ice  by  meas¬ 
urements  at  different  frequencies.  We  pick  an  overly 
simplified  case  for  illustration. 

Let  us  suppose  that  the  various  components  of  a  sam¬ 
ple,  ice,  water  and  air,  could  be  modeled  as  acting  in 
parallel  as  shown  in  Figure  3,  a  unit  cube.  If  electrodes 
are  applied  at  the  top  and  the  bottom,  the  contributions 
to  the  measured  capacitance  of  each  material  would  be 
in  parallel.  Let/A,/w.  and/j,  be  the  volume  fractions  of 
each  material  and  C^,  and  C\  the  partial  capacitances 
due  to  air,  water  and  ice.  respectively.  Using  eg  as  the 
permittivity  of  free  space  (air)  and  and  k\  as  the 
dielectric  constants  of  water  and  ice,  we  have  C  mea¬ 
sured  (Cm)  given  by 


Cm  —  Ca  +  Cw -!■  Cj  —  Eo(/a  +  ^wi/w  ^j^i)-  (1) 

Also 

/a  +/w  +/a  =  1  (unit  cube).  (2) 

Since  k\  (high  frequency)  is  very  different  from  k\ 
(low  frequency),  if  we  measure  C  at  both  high  and  low 
frequencies,  we  have  three  different  independent  equa¬ 
tions  and  three  unknowns.  Therefore,  we  should  be  able 
to  solve  for  /a,  /w  and  /j,  without  an  independent 
measurement  of  density.  Unfortunately  the  situation  is 
more  complicated  for  two  basic  reasons.  The  first  is 
electrical  and  has  already  been  touched  on.  The  second 
is  geometrical. 

The  electrical  problem  with  capacitance  measure¬ 
ment  is  that  the  low  frequency  capacitance  is  complicat¬ 
ed  by  electrode  and  impurity  effects  and  the  possibility 
of  very  low-frequency  dispersions  for  ice,  which  have 
been  reported  in  the  literature  (Camp  et  al.  1967,  Von 


Hippel  1972).  Perhaps  this  can  be  overcome  by  making 
measurements  at  several  different  frequencies  through¬ 
out  the  Debye  dispersion  range  and  using  these  to  effec¬ 
tively  project  the  true  low-frequency  value. 

To  illustrate  the  geometrical  problems,  we  refer 
again  to  Figure  3.  If  instead  of  placing  electrodes  on  the 
top  and  bottom  of  the  cube,  we  place  them  on  the  left- 
hand  and  right-hand  sides,  we  have  an  entirely  different 
geometry.  The  air,  water  and  ice  are  now  in  series  in¬ 
stead  of  in  parallel.  Our  measurements  would  give 
entirely  different  values,  as 

i/CA+i/c^+i/q 


instead  of  Ca  +  +  Cj  in  the  parallel  case. 

As  a  specific  illustration,  let  us  suppose  /a  =  0.8,/;  = 
0. 1 5  and/,v  =  0.05.  For  the  parallel  case  (eq  1 ),  using 
(low)  =  100,  k^  =  88,  and  k^  (high)  =  3.5, 

Cmeas  (low  freq)  =  20.2eo 

Cmeas  (high  freq)  =  5.73eo. 


In  the  series  case  ,  we  would  have, 


^^rneas  (series)  =£0 


/A  +  /i/^i+/«/^v 


(3) 


Cmeas  (low  freq)  =  1.247eo 

Cmeas  (high  freq)  =  1.186eo. 

These  approximate  the  two  extreme  cases  for  the  possi¬ 
ble  grain  structure  of  snow,  parallel  columns  and  series 
plates.  Thus,  from  the  measurements  themselves,  one 
might  be  able  to  infer  something  about  the  grain  struc¬ 
ture.  However,  the  situation  is  even  more  complicated 
than  this.  Temperate  snow  is  a  complex  stmeture  of 
grains  of  ice  of  various  shapes  in  series  and  parallel 
combinations  with  water  that  has  electrical  conductiv¬ 
ity  as  well  as  dielectric  constant.  (This  is  also  a  problem 
with  microwave  measurements.)  Thus  in  the  end,  the 
measurements  must  be  calibrated  against  samples  hav¬ 
ing  known  properties.  But  these  complicating  factors 
have  their  positive  side  as  well.  If  we  measure  at  a  num¬ 
ber  of  points  over  the  range  of  the  Debye  dispersion,  we 
might  be  able  to  develop  a  set  of  profiles  each  of  which 
is  uniquely  characteristic  of  a  specific  water  content, 
density  and  snow  structure.  The  studies  reported  here 
were  undertaken  to  explore  this  possibility. 

In  this  specific  effort,  we  have  been  so  far  unsuccess¬ 
ful.  However,  in  the  process,  we  have  developed  what 
we  believe  is  a  significant  tool  for  studies  of  wet  snow. 
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some  useful  insights  into  the  process  involved,  and  a 
body  of  specific  information  concerning  the  dielectric 
properties  of  snow,  both  wet  and  dry,  in  this  frequency 
region.  In  order  to  make  this  information  generally 
available  and  useful  we  give  in  this  report  a  rather  de¬ 
tailed  account  of  our  experiments  and  their  results.  We 
have  found  no  such  extensive  information  in  the  liter¬ 
ature  and  hope  that  this  work  will  provide  a  useful  basis 
for  theoretical  or  further  experimental  studies. 

BACKGROUND 
Literature  review 

A  computer  search  was  conducted  covering  the  fol¬ 
lowing  databases  and  years:  “Information  Service  for 
Physics  and  Engineering  (INSPEC),”  1969-1989;  and 
“Engineering  Index,”  1970-1989.  This  was  supple¬ 
mented  by  a  manual  search  of  the  CRREL  abstracts  for 
the  years  available  to  us,  1961-75.  We  are  indebted  to 
Ronald  Atkins  who  has  provided  us  with  further  titles 
from  the  CRREL  abstracts.  All  references  that  we  have 
found  useful  to  this  line  of  study  are  included  in  the 
Annotated  Bibliography. 


Modeling  by  computer 

Because  any  electrical  analog  we  might  make  in¬ 
volves  the  series-parallel  interconnection  of  resistors 
and  capacitors,  an  interactive  computer  program  has 
been  developed  (using  Microsoft  Quick-Basic,  version 
4)  that  allows  us  to  model  resistance  to  capacitance 
circuits  of  almost  any  complexity.  If  we  have  a  model  in 
mind,  we  can  try  it  out  and  plot  its  capacitance  and  con¬ 
ductance  spectrum.  We  can  also  go  to  the  other  extreme 
and  find  out  how  tinkering  with  a  circuit  will  change  its 
characteristics.  Figures  4  and  5  are  examples  of  the  first 
process.  We  happened  to  have  available  dielectric  spec¬ 
tra  for  Portland  cement  samples  and  a  model  that  had 
been  proposed  to  describe  two  different  states  of  curing. 
By  adjusting  parameters,  we  were  able  to  fit  the  data  for 
cement  cured  only  one  day  (Fig.  4)  and  cement  cured  35 
days  (Fig.  5).  Figures  6  and  7a  are  examples  of  the 
second  process.  We  asked  whether  it  was  possible  to 
create  from  resistors  and  capacitors  alone  a  network  that 
would  have  a  dispersion  curve  in  which  the  capacitance 
would  drop  (and  the  conductance  would  rise)  over  a 
wide  frequency  range,  not  as  frequency  squared,  but  as 
frequency  to  a  lower  power  (as  happens  in  snow).  These 
figures  show  that  it  can  be  done.  Figure  6  shows  that  by 


Log  Cp(F)  Log  Gp(S) 


Figure  4.  E.xample  of  matching  a  model  to  experimental  results  using  data  available  from  dielectric 
spectra  of  portland  cement  paste  cured  for  one  day. 
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Figure  5.  Example  of  matching  a  model  to  experimental  results  using  data  available  from  dielectric 
spectra  of  portland  cement  paste  cured  for  35  days. 


Figure  6.  Equivalent  parallel  capacitance  and  conductance  as 
functions  of frequency  for  the  network  shown. 


5 


Log  Cp(F)  Log  Gp(S) 


a.  Equivalent  parallel  capacitance  and  conductance  as  functions 
of  frequency  for  the  nework  shown. 


K' 

b.  Cole  plot  for  the  network  of  7a. 
Figure  7.  Properties  of  the  network  shown. 


doubling  C  per  section  and  keeping  G  coi .slant,  the  total 
capacitance  can  be  made  to  fall  as  the  first  power  of  the 
frequency  from  about  200  Hz  to  20  kHz.  Figure  7a 
shows  that  if  C  is  halved  and  G  is  doubled  for  each 
succeeding  section  the  total  C  falls  and  the  total  C  rises 
as  the  square  root  of  frequency  over  about  three  de- 
cades.  Figure  7b  shows  how  the  Cole  semicircle  is  dis¬ 


torted  for  this  distribution  of  Debye  relaxations.  Both  of 
these  features  show  up  in  our  data  for  snow. 

Computer  tools  for  data  management 
and  interpretation 

Each  snow  experiment  involved  ihe  collection,  pro¬ 
cessing  and  organization  of  a  large  variety  of  data. 
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These  data  included  time,  sample  temperature,  plate 
separation,  snow  parameters,  power  measurement  and 
dielectric  properties  as  a  function  of  frequency.  The 
data  entry,  processing  and  organization  was  simplified 
by  the  use  of  a  computer  spreadsheet  (Quattro  by  Bor¬ 
land)  on  a  standard  IBM  AT  microcomputer.  The  spread¬ 
sheet  was  used  for  the  calculation  of  snow  water  con¬ 
tent.  During  the  dielectric  heating  of  the  snow,  the 
power  was  calculated  from  the  current,  voltage  and 
phase  measurements.  Numerical  integration  of  power 
during  the  heating  process  determined  the  heat  input 
and  thus  the  water  content.  The  spreadsheet  was  also 
used  in  calculating  the  conductivity  and  dielectric  con¬ 
stants,  and  the  data  for  the  logarithmic  graphs  and  the 
Cole  plots. 

To  further  analyze  the  data,  an  interactive,  curve¬ 
fitting  program  was  written  in  Microsoft  Quick-Basic, 
version  4.This  program  finds  the  parameters  giving  the 
best  fit  for  a  set  of  data  to  any  type  of  polynomial  speci¬ 
fied  by  the  user. 

EXPERIMENTAL  DESIGN 

Dielectric  measurements 

Measurements  were  made  of  equivalent  parallel 
capacitance  Cp  and  equivalent  parallel  conductance  Gp 
using  a  General  Radio  1615A  bridge  in  the  three- 
terminal  mode.  The  sample  holders  are  shown  in  Fig¬ 
ures  8  and  9.  The  signal  source  was  a  General  Radio 
1210-C  oscillator  and  the  detector  was  the  General 
Radio  1232- A  tuned  amplifier  and  null  detector.  The 
sample  cells  were  placed  inside  an  electrically  shielded 
oven,  which  was  itself  inside  a  commercial  chest  type 
freezer.  The  oven  temperature  was  controlled  by  a 
Cole-Palmer  model  2158  Versatherm  controller  using 
a  thermistor  sensor.  Sample  and  oven  temperatures 
were  monitored  by  Yellow  Springs  Instrument  linear¬ 
ized  thermistors  with  digital  readout  on  Keithley  model 
160  voltmeters.  All  probes  were  calibrated  against  a 
Hewlett-Packard  2802-A  platinum  resistance  thermo¬ 
meter  and  zeroed  in  a  common  ice  bath. 

Sample  water  content 

An  obvious  way  to  prepare  snow  of  known  water 
content  is  simply  to  add  to  dry  snow  at  0°C  known 
amounts  of  water  also  at  0°C.  However,  there  are  prob¬ 
lems  with  this  method.  One  of  these  is  obtaining  uni¬ 
form  mixing.  Another  is  that  if  either  the  snow  or  the 
water  is  too  cold,  some  of  the  water  will  freeze.  If  the 
water  is  slightly  too  warm,  some  of  the  ice  will  melt. 
Most  importantly  each  time  the  water  content  is  changed, 
we  must  work  with  a  new  and  possibly  fundamentally 
different  sample  that  inevitably  will  be  mixed  and 
packed  differently.  What  we  would  like  to  do  is  to  vary 


the  water  content  of  a  particular  sample  which  remains 
undisturbed  in  the  test  cell.  To  this  end,  we  have  devel¬ 
oped  a  method  of  electrically  melting  the  sample  itself 

In  doing  this,  two  things  are  important;  the  heat  gen¬ 
erated  must  be  uniformly  deposited  in  the  sample  and  it 
must  be  measurable  with  sufficient  precision.  Direct 
current  heating,  perhaps  the  simplest  approach,  is  unsat¬ 
isfactory  because  once  a  liquid  path  has  been  created, 
ionic  conduction  will  cause  most  of  the  heat  to  be  de¬ 
posited  along  that  path,  which  will  become  larger  and 
larger.  To  ensure  that  the  energy  is  uniformly  deposited 
in  the  ice  itself,  we  use  the  heating  due  to  dielectric  loss 
in  the  ice.  Thus  we  choose  a  frequency  of  the  order  of, 
or  higher  than,  the  Debye  dispersion  frequency  for  ice 
(/d).  If  we  used  microwave  heating,  we  would  again 
find  water  absorbing  much  of  the  energy.  Moreover, 
microwave  heating  would  be  hard  to  monitor. 

We  have  chosen  20  kHz  as  our  working  frequency. 
The  power  is  monitored  by  measuring  the  voltage  ap¬ 
plied  to  the  sample,  the  current  through  it  and  the  phase 
angle  between  them.  A  Singer  model  LTPS-Il  oscillator 
and  power  amplifier  supplies  up  to  1600  V  at  20  kHz 
from  a  780-ohm  source.  Thus  arbitrary  increments  of 
heat  can  be  added  to  a  sample  in  si  tu  and  dielectric  meas¬ 
urements  can  be  made  after  each  increment  without 
disturbing  the  sample.  The  circuit  arrangement  is  shown 
in  Figure  10. 

The  most  difficult  aspect  of  this  procedure  is  the 
measurement  of  the  phase  angle,  which  changes  signifi¬ 
cantly  with  time.  At  the  start,  we  monitored  this  using  a 
dual  trace  oscilloscope.  More  recently,  we  have  been 
provided  by  CRREL  with  a  Wavetek  model  750  phase 
meter.  We  now  use  both.  In  two  cases,  we  checked  our 
water  content  as  calculated  by  power  measurement 
against  alcohol  calorimetry.  In  one,  the  energy  calcula¬ 
tion  gave  a  20%  water  content  and  the  alcohol  calori¬ 
meter  a  value  of  23%.  In  the  other  the  energy  calculation 
gave  5%  and  the  alcohol  calorimetry  7.8%.  To  check 
further  on  the  alcohol  calorimetric  technique,  we  sim¬ 
ply  added  5%  and  1 0%  water  to  snow  samples  at  0°C 
and  found  that  alcohol  calorimetry  gave  us  7.4%  and 
14.5%  readings  respectively.  This  error  was  consistent 
with  the  higher  readings  given  above  and  gives  us 
greater  confidence  in  our  power  measurement  proce¬ 
dure  than  in  our  alcohol  calorimetry. 

We  appreciate  the  comment  of  a  reviewer  who  has 
called  our  attention  to  the  fact  that  E.  Brun  (1989)  had 
previously  used  dielectric  heating  at  20  kHz  to  change 
the  water  content  of  snow  by  melting.  In  that  paper, 
Brun  did  not  discuss  his  measurement  of  power. 

Monitoring  dimensional  change 

In  preliminary  experiments,  we  used  a  capacitance 
cell  that  had  vertical  plates  1 0  cm  square  separated  by  2 
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Figure  8.  Measurement  cell  for  experiments  I,  2  and  3.  The  electrodes  are  thin 
copper  foil. 


cm.  A  vertical  structure  was  used  so  that  if  water  drained 
away  from  the  top  of  the  sample,  good  electrical  contact 
would  be  maintained  over  most  of  the  sample  and  a 
good  average  measurement  would  still  be  made.  How¬ 
ever,  we  found  that  even  a  very  small  amount  of  mois¬ 
ture  caused  the  sample  to  shrink  away  from  the  plates. 
This  we  attribute  to  capillarity  at  low  water  contents.  At 
high  water  contents,  there  was  also  shrinkage  due  to  the 
melting  of  ice.  Thus  we  changed  our  basic  design  to  a 
horizontal  cell  (Fig.  8)  and  included  in  it  a  device  for 
monitoring  the  sample  thickness  by  mechanically  link¬ 
ing  the  top  plate  to  a  high  resolution  potentiometer  on 
top  of  the  sample  cell.  The  resistance  of  the  potentiome¬ 
ter  was  monitored  on  a  chart  recorder,  which  provided 


a  continuous  record  of  the  position  of  the  top  plate.  This 
also  allowed  us  to  correct  for  changes  in  sample  thick¬ 
ness  when  calculating  conductivity  and  dielectric  con¬ 
stant. 

Overall  design 

Except  where  noted,  all  the  measurements  reported 
here  were  made  using  one  of  two  sample  cells.  The  first 
is  shown  in  Figure  8.  The  housing  is  made  of  building- 
grade  dense  Styrofoam  having  an  R  value  of  6  for  a  1- 
in.  thickness  (thermal  conductivity  of  0.024  W/m  K). 
The  electrodes  are  of  thin  copper  foil  to  minimize  the 
heat  capacity.  They  are  unguarded  for  flinging  cap¬ 
acitance.  The  guard  terminal  of  the  bridge  is  connected 
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Figure  9.  Measurement  cell  for  experiments  10-16.  In  experiment  10,  segments  2, 3  and 
4  were  guards  and  in  experiments  11-16,  segments  3  and  4  were  guards  and  1  and  2 
connected  together. 


to  the  aluminum  case  of  the  oven.  On  top  of  the  upper  temperature  of  the  sample.  Figure  9  shows  the  design  of 

electrode  and  insulated  from  it  both  thermally  and  elec-  the  second  cell  that  we  used.  This  was  a  guarded  cell 

trically  is  a  brass  disk  that  provides  a  constant  gravita-  using  V4-in.  stainless  steel  electrodes.  The  loading  pro- 

tional  loading  of  1027  N/m^  (about  0. 1 5  Ib/in.^).  This  is  vided  by  this  cell  was  306  N/m^.  In  all  other  respects  it 

done  to  ensure  that  all  samples  in  these  experiments  was  very  similar.  (We  show  later  that  loading  itself  has 

have  the  same  packing,  thus  eliminating  one  possible  a  negligible  effect  on  our  measurements  except  as  it 

variable.  A  rod  from  this  plate  drives  the  precision  changes  density.  The  density  of  each  sample  was  sepa- 

potentiometer  to  provide  displacement  measurements  rately  measured.)  Each  cell  itself  was  placed  inside  a 
as  the  sample  thickness  changes.  A  thermistor  bead  in  thick  walled  aluminum  “oven,”  which  was  in  turn  inside 

contact  with  the  bottom  of  the  bottom  plate  monitors  the  the  commercial  freezer.  The  oven  was  wound  with  re- 
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Position  1  Heat 

2  Measure 


Figure  10.  Circuit  arrangement  for  measuring  dielectric  properties  of  wet  snow. 


sistance  wire  through  which  DC  current  was  passed  to 
maintain  the  desired  temperature.  This  was  controlled 
by  a  thermistor  sensor  and  Cole-Palmer  Model  2158 
temperature  controller.  The  oven  temperature  was  inde¬ 
pendently  monitored  and  was  normally  maintained  at 
-1  (±0.5)  °C. 

Figure  10  is  a  diagram  of  the  electrical  connections 
to  the  test  cell.  Stray  capacitances  due  to  shielding  of  the 
leads,  etc.,  are  shown  in  dotted  lines.  Since  the  source  is 
of  low  impedance,  C I  has  negligible  effect.  Because  of 
the  bridge  design,  C2  and  C3  can  be  neglected.  Because 
/?!  is  so  small,  it  has  a  negligible  effect  on  the  measure¬ 
ment  of  y  and  C4  of  this  magnitude  can  be  neglected. 

The  performance  characteristics  are  as  follows:  Cf 
and  Gf  are  ordinarily  measured  to  4  or  more  significant 
figures  with  an  accuracy  of  0. 1  %  over  most  of  the  range. 
However,  to  translate  this  into  dielectric  constant  and 


conductivity  requires  an  accurate  knowledge  of  sample 
dimensions,  which  are  used  in  calculating  the  cell 
capacitance  and  the  effects  of  fringing.  The  dominant 
uncertainty  is  in  sample  thickness.  We  estimate  a  com¬ 
bined  uncertainty  of  10%  in  absolute  value  and  a  reso¬ 
lution  of  better  than  1  %.  Density  is  measured  by  weigh¬ 
ing  the  sample  to  better  than  1%  and  measuring  its 
thickness  after  the  cell  has  been  assembled.  The  abso¬ 
lute  accuracy  is  estimated  to  be  10%,  the  accuracy  of 
changes  about  0. 1  %.  We  can  resolve  a  change  of  thick¬ 
ness  of  0.01  mm  while  thickness  is  accurate  to  1  mm. 
Temperatures  are  ordinarily  measured  with  an  accuracy 
of0.2®Cand  a  resolution  ofO.OJ  °C.  Overall  we  believe 
our  water  content  is  accurate  to  0. 1  times  the  reading 
±0.5%,  and  our  values  for  conductivity  and  dielectric 
constant  are  accurate  to  10%  of  the  reading  with  a 
resolution  of  better  than  1  %. 


Table  1.  Properties  of  snows  used. 


Harvest  data _  Melt _ Analv.sis  In)  (eg  I) 


Snow 

Date 

Temp 

(gicnr') 

pH 

0* 

Ca 

Mg 

K 

Na 

NH4 

a 

NOj 

SO4 

lA 

March  89 
granular 

0.42 

4.5 

169 

IB 

March  89 
Granular 

0.52 

6.04 

3.56 

227 

10.7 

4.3 

74.4 

1.6 

70 

10.5 

203 

11 

30  Dec  89 
very  fluffy 

-8 

0.38 

4.27 

2.80 

11.5 

3.3 

4.3 

52.6 

nd 

56 

45.8 

21 

III 

30  Jan  90 

_2 

0.41 

5.36 

0.305 

5.0 

nd*' 

1.5 

1.7 

nd 

3 

19 

5 

IV 

5  Feb  90 
light  fluff 

-15.5 

0.32 

5.0 

0.534 

7.0 

nd 

1.3 

1.3 

nd 

4 

8.2 

6 

V 

25  Feb  90 
fine  &  floe. 

-8 

0.27 

4.82 

0.865 

19.0 

1.6 

4.3 

7.4 

nd 

II 

13.4 

17 

VI** 

31  Mar  90 

0 

5.64 

0.17 

0.1 

nd 

4.8 

nd 

nd 

6 

0.6 

3 

wet 


*  Density  as  measured  on  14  April  1990. 
t  Units  of  conductivity  arc  milli  (ohm-meters)  -I  or  milli  S/m. 
*♦  Snow  VI  not  used  in  measurements  reported  here, 
tt  No  data. 
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S  AMPl.i;  l»Ri;i‘  AR  A  HON  AM)  I’RORKR  riKS 
Harxcsf  and  stnrajic 

Sn(".\  samples  uerecDlleeled  loeall>  and.  e\eepi  tor 
snow  I.  snoii  allei  llie  storm  liad  einied.  .\t  llte  time  ol 
eolleeiion,  tlie  sno\\  iem|K'raUne  v\  as  measured  and  the 
sirueliire  noted,  llte  snou  utis  jilaceti  in  .Si\roloam 
l'o\es  and  storerl  m  a  tiee/er  at  I  2  ('  unlil  it  w  ;is  used. 
Tiie  properties  ot  tlie  snow  s  used  are  listed  in  Table  I 
Miero|iliotoi:taphs  ol  snows  I  ihroueli  \  taken  on  S 
Mareli  !')')()  are  shown  in  l  ieure  1  1.  Snow  I  lA  ami  Bi 
was  haiwested  in  Mareli  IdXh.  It  had  been  on  the  eround 
lor  more  than  a  week  aiul  alreade  IkkI  undergone  eon 
siderable  metamorphosis  when  h.n'\  ested.  It  w  as  oriui- 
nalB  stored  in  poK ethylene  hues  in  the  Iree/er  anil  w  as 
Inst  iisetl  m  September  I  bSM.  The  others  were  hare  ested 
in  late  IdS')  ami  earh  Iddd  as  mdieated  in  Table  1 


Properties 

One  or  more  samples  of  each  snow  harvest  were 
melted  and  the  conductivity  and  pH  of  each  determined. 
In  addition  an  impurity  analysis  was  made.  The  results 
are  summarized  in  Table  1.  Not  surprisingly,  there  is  a 
strong  correlation  between  the  NaCl  concentration  and 
the  conductivity  of  the  melt.  It  may  also  be  of  interest 
that  the  charge  of  the  total  anions  minus  the  total  cations 
is  positive  for  the  two  lowest  pH  snows  and  negative  for 
that  having  the  highest  value. 

Since  the  same  snow  harvest  was  used  for  a  number 
of  experiments  at  different  times  and  since  properties 
such  as  density  and  grain  size  vary  with  time,  each 
specific  sample  is  identified  by  an  (arabic)  experiment 
number.  Table  2  is  an  index  to  the  samples. 

Table  2.  Sample  data. 


Measurement 


Experiment 

date 

Snow 

Code 

Experiment 

1 

5  Sept  89 

lA 

8/89-A 

K<5%W 

2 

1 1  Sept  89 

lA 

9/89-B 

Ka%W 

3 

15  Sept  89 

lA 

8/89-C 

Ka%W 

I/90-A 

Test  run  of 

no  interest 

1/90-B 

Test  run  of 

no  interest 

4 

24  Jan  90 

II 

I/90-C 

Density 

5 

31  Jan  90 

III 

1/90-D 

Density 

Loose  guard  ring 

6 

31  Jan  90 

III 

1/90-E 

Density 

7 

5  Feb  90 

tv 

2/90-A 

Density 

8 

8  Feb  90 

II 

2/90-B 

Density 

9 

12  Feb  90 

IB 

2/90-C 

Density 

10 

1  Mar  90 

V 

3/90-A 

K<s%W 

II 

6  Mar  90 

IV 

3/90-B 

K<s%W 

12 

7  Mar  90 

II 

3/90-C 

K<s%W 

13 

8  Mar  90 

IB 

3/90-D 

Ko%W 

14 

9  Mar  90 

III 

3/90-E 

K(s%W 

15 

15  Mar  90 

II 

3/90-F 

Krs%W 

16 

20  Mar  90 

IB 

3/90-G 

Ko%W 

CONTROL  EXPERIMENTS 

At  various  times  during  our  experiments,  questions 
arose  about  possible  sources  of  error.  Each  time  we  tried 
to  devise  experiments  to  check  these.  Although  they 
were  conducted  as  the  questions  arose,  they  are  grouped 
together  in  this  section. 

Electrode  materials 

A  recurrent  question  in  our  early  experiments  was 
whether  part  of  the  dispersion  we  were  seeing  was  an 
electrode  effect  characteristic  of  the  electrode  material. 
To  explore  this,  we  varied  electrode  materials.  A  prob¬ 
lem  with  Just  using  cells  of  various  materials  is  that  then 


we  would  always  be  comparing  different  samples.  Thus 
we  used  a  cell  with  a  bottom  plate  of  copper  and  changed 
the  top  plate  without  disturbing  the  sample.  Top  plates 
of  aluminum,  copper,  stainless  steel  and  a  blocking 
electrode  formed  by  insulating  the  top  plate  with  a  very 
thin  layer  of  polyethylene  were  used.  Measurements 
were  made  from  50  Hz  to  50  kHz  with  very  similar 
results  in  each  case.  The  maximum  difference  in  Cp, 
which  was  found  at  50  Hz,  was  only  12%.  That  in  Gp 
was  only  8%.  We  conclude  that  the  electrode  contribu¬ 
tion  was  not  a  serious  problem. 

Sample  distribution  in  the  cell 

Another  concern  was  that  our  results  might  be  signif¬ 
icantly  influenced  by  an  undetected  variation  of  the 
distribution  of  snow  in  the  cell.  Snow  was  typically 
placed  in  the  cell  as  evenly  as  possible  and  then  leveled 
with  a  blade.  To  test  the  distribution  we  built  the  cell  of 
Figure  9  using  stainless  steel  electrodes  (this  became 
our  cell  of  choice  for  later  experiments).  The  top  elec¬ 
trode  was  one  piece,  10  cm  in  diameter.  The  bottom 
plate  comprised  a  central  disk  surrounded  with  three 
concentric  rings  with  about  one  millimeter  of  insulation 
between  them.  The  central  disk  and  each  ring  had  an 
area  of  20  cm^.  By  connecting  all  of  the  elements  but 
one  to  the  guard  circuit,  we  were  able  to  measure  each 
region  of  the  sample  in  turn  without  disturbance. 

Several  experiments  with  snow  and  one  using  wet 
sand  suggested  a  possible  small  radial  decrease  in  Cp 
and  Gp,  with  the  largest  variation  at  low  frequencies. 
This  induced  us  to  use  greater  care  in  uniformly  filling 
the  cells. 

We  also  considered  the  possibility  of  an  azimuthal 
variation  of  snow  distribution.  To  test  for  this  possibil¬ 
ity,  we  made  a  top  plate  that  was  segmented  like  a  pie. 
We  were  able  to  sample  one  segment  at  a  time  and  found 
no  significant  azimuthal  variation.  This  test  also  ver¬ 
ified  the  uniformity  of  the  test  samples  themselves. 

Effects  of  pressure 

As  explained  earlier,  a  standard  sample  loading  was 
used  to  eliminate  possible  variations  due  to  pressure 
effects.  However,  the  following  experiment  was  done  to 
explore  the  matter.  A  sample  was  measured  with  the 
standard  load.  Then  incremental  loads  of  5(X)  g  were 
added  and  the  sample  remeasured  with  each  addition. 

Each  time  the  samples  became  slightly  more  com¬ 
pacted  and  the  values  of  Cp  and  Gp  increased  by  2-3%. 
But  when  all  the  extra  loading  was  removed,  the  values 
remained  almost  unchanged.  We  conclude  that  the 
observed  increases  were  due  to  increases  in  density,  not 
pressure.  Later  experiments  on  the  effect  of  density  (VI 
B)  confirm  the  large  effect  of  density. 


Figure  12.  Temperature,  power  le\’el,  and  sample  thickness  as  a  function  of  time 
during  a  typical  experiment 


DIELECTRIC  EXPERIMENTS  FOR 
WET  AND  DRY  SNOW 

Experiments  on  aged  snow 

Because  it  was  late  December  1989  before  we  were 
able  to  harvest  fresh  snow,  many  of  our  experiments  had 
to  use  snow  that  had  been  harvested  in  the  previous 
March  (snow  I).  This  snow  had  already  been  on  the 
ground  for  a  week  or  more  before  it  was  collected.  Con¬ 
siderable  metamorphosis  had  already  taken  place  when 
it  was  first  used.  Its  characteristics  are  listed  in  Table  1 . 
Figure  1 1  is  a  microphotograph  of  its  structure  taken  in 
March  of  1 990  when  it  was  almost  a  year  old.  At  the  time 
of  these  measurements,  a  typical  grain  dimension  was 
about  2  mm.  For  these  experiments,  cell  I  (Fig.  8 )  was 
used. 

Figure  1 2,  a  plot  of  temperature,  power  and  sample 
thickness,  illustrates  how  the  experiments  were  con¬ 
ducted.  The  sample  was  placed  in  the  cell  at  about  -6°C, 
which  was  placed  in  the  “oven”  at  -0.5°C.  At  this  time, 
a  measurement  of  Cp  and  Gp  as  a  function  of  frequency 
(a  frequency  scan)  was  made. 

Over  the  first  32  minutes,  the  sample  warmed  slowly 
to  about  -3°C  at  which  point  (point  A),  the  20-kHz 
power  was  turned  on  and  left  on  until  B,  heating  the 
sample  internally  almost  to  zero.  A  second  frequency 
scan  was  then  made.  At  about  62  minutes  (point  Q,  the 
heating  power  was  again  turned  on  briefly  until  D. 
During  this  time  ice  was  melted.  Again  a  frequency  scan 
was  run,  followed  by  another  heat  pulse,  etc.  Note  the 
gradual  rise  in  the  sample  thickness  for  about  the  first  30 


minutes  as  if  the  sample  were  relaxing  slightly  after  be¬ 
ing  loaded.  Then  the  gradual  packing  at  zero  degrees, 
followed  by  a  sharp  decrease  at  C  as  ice  melts  followed 
by  successive  decreases  at  E  and  G.  These  last  two  de¬ 
creases  do  not  stop  abruptly  with  thecessation  of  power. 
Moreover,  the  rate  of  drop  is  not  as  great  as  the  first, 
even  though  the  power  is  larger.  We  suggest  that  capil¬ 
lary  forces  may  play  a  significant  role  in  the  first  and  that 
the  others  are  dominated  by  volume  loss  due  to  melting 
and  rearrangement  of  the  grains.  The  situation  is  illus¬ 
trated  in  Figure  1 3  in  which  the  electrode  separation  is 
shown  to  change  as  the  water  content  is  changed  by 
heating. 

Data  for  samples  1,2  and  3  are  shown  in  Figures  14- 
1 6.  These  measurements  were  all  made  using  cell  I  with 


Figure  13.  Electrode  separation  vs  percent  water 
for  sample  3,  which  shows  the  sample  shrinking 
as  percent  water  increases. 
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iBX. 


copper  foil  electrodes.  The  curves  for  the  log  lo  conduc¬ 
tivity  vs  log  10  frequency  change  their  curvature  in  going 
from  dry  snow  to  even  slightly  damp  snow.  Then  the 
predominant  change  is  a  steady  rise  of  the  values  with 
water  content.  The  increase  with  frequency  is  slow  and 
almost  linear.  The  curves  for  the  logio  of  apparent  di¬ 
electric  constant  k  also  show  a  sudden  change  in  going 
from  dry  snow  to  wet  snow.  Here  again  there  is  strong 
sample-to-sample  similarity  except  for  the  highest  fre¬ 
quency  points  of  sample  1 .  The  log  lo  of  falls  roughly 
in  proportion  to  the  logiQof  frequency.  It  was  shown  in 
the  Modeling  by  Computer  section  that  this  type  of  dis¬ 
persion  in  both  the  dielectric  constant  and  conductivity 
is  characteristic  of  a  system  having  a  number  of  Debye 
relaxations  in  parallel,  each  having  a  different  relaxa¬ 
tion  time.  This  is  also  seen  in  the  Cole  plots  (Fig.  14- 
16c-f). 

We  have  used  above  the  term  “apparent  dielectric 
constant”  because  at  low  frequencies  we  frequently  find 
values  greatly  exceeding  that  for  pure  ice  or  pure  water. 
Thus  the  microscopic  geometry  of  the  sample  must  be 
very  different  from  the  geometry  of  our  cell.  There  may 
be  conducting  paths  that  make  the  effective  plate  area 
much  greater  than  the  tme  plate  area  and  the  effective 
plate  separation  small.  By  the  same  argument,  the  term 
“effective  conductivity”  might  be  used.  In  the  interest  of 
brevity,  we  will  drop  the  word  “effective”  in  the  re¬ 
mainder  of  this  report. 

The  data  for  these  three  samples,  all  from  snow  lA, 
form  a  consistent  set  for  10  different  water  contents.  For 
zero  water  content,  they  result  in  almost  identical  spec¬ 
tra,  giving  us  confidence  in  the  reproducibility  of  our 
experiments.  When  the  water  contents  of  different  ex¬ 
periments  are  close,  the  spectra  are  very  similar.  This 
tends  to  confirm  our  estimates  of  the  reliability  of  the 
water  content  (±  1  %).  They  were  all  taken  within  a  space 
of  1 0  days.  The  data  for  samples  1 3  and  1 6  (Fig.  1 7  and 
18)  are  for  snow  sample  IB,  which  was  harvested  with 
sample  lA  but  from  a  different  region  and  has  a  different 
conductivity  and  pH.  These  data  were  taken  six  months 
later  than  those  for  snow  lA  and  with  the  cell  of  Figure 
9.  The  figures  are  very  different  indeed,  a  fact  that  is 
discussed  in  the  Discussion  below.. 

Experiments  on  new  snow 

After  the  winter  snowfall  began,  we  were  able  to  ex¬ 
periment  with  different  snows.  For  all  of  these  we  used 
the  stainless  steel  cell  of  Figure  9.  The  center  plate  and 
first  ring  were  connected  together  as  one  electrode  and 
the  outer  two  electrodes  connected  to  the  guard.  (In  the 
case  of  sample  10,  all  three  outer  rings  were  connected 
to  guard.)  The  procedure  was  essentially  the  same  as 
that  used  in  A.  First,  frequency  runs  were  made  on  the 
sample  when  it  was  several  degrees  below  zero,  then  the 


sample  was  heated  to  very  close  to  zero  using  20  kHz,  as 
explained,  and  measured  again.  Thereafter,  increments 
of  heat  were  added  to  melt  some  snow,  and  frequency 
scans  were  run  after  each  increment.  Figures  17-23 
show  the  data  for  these  samples  both  as  plots  of  the  log 
conductivity  and  logio  dielectric  constant  vs  log  fre¬ 
quency  and  as  Cole  diagrams.  These  data,  while  cover¬ 
ing  the  same  general  range  of  values,  differ  in  detail 
from  one  another  and  from  the  earlier  measurements. 

Discussion 

Our  major  experimental  results  are  summarized  in 
Figures  14-23.  Numerical  values  are  tabulated  in  Ap¬ 
pendix  A.  For  detailed  information  concerning  the  prop¬ 
erties  of  any  particular  sample,  see  Tables  1  and  2.  The 
dominant  source  of  absolute  error  in  the  measurement  of 
conductivity  and  dielectric  constant  is  the  sample  thick¬ 
ness.  This  is  less  than  10%.  However,  the  conductance 
and  capacitance  were  measured  to  0.1%  or  better.  We 
believe  that  the  major  source  of  error  in  the  water  content 
was  inaccuracy  in  the  power  measurement.  This  might 
cause  an  error  between  samples  such  that  one  sample 
measured  at  8%  might  have  the  same  water  content  as 
another  measured  at  10%.  But  measurements  on  the 
same  sample  should  be  ordered  so  that  the  relative  errors 
are  considerably  less. 

The  general  features  of  the  dielectric  data  for  all 
samples  show  a  rather  large  low-frequency  conductivity 
rising monotonically  with  frequency.  This  is  compatible 
with  a  conductivity  made  up  of  two  parts,  a  frequency- 
independent  ionic  conductivity  (protons)  to  which  is 
added  a  conductivity  resulting  from  one  or  more  dielec¬ 
tric  dispersions.  The  dielectric  constant  data  show  these 
dispersions  more  clearly  and  indicate  a  range  of  disper¬ 
sions.  The  Cole  plots  show  a  range  of  dispersions  in  the 
low-frequency  range  centered  in  the  region  of  100-200 
Hz  and  often  a  range  of  dispersions  in  the  kilohertz 
range.  In  general,  both  the  conductivity  and  the  dielec¬ 
tric  constant  rise  with  water  content.  On  many  samples, 
there  is  a  change  in  curvature  from  concave  up  to  con¬ 
cave  down  on  the  log  jq  conductivity  and  log  iq  dielectric 
constant  vs  frequency  plots  as  the  sample  changes  from 
dry  to  wet. 

For  dry  snow,  there  appears  to  be  a  general  tendency 
for  the  low-frequency  conductivity  to  be  proportional  to 
the  H'*'  ion  concentration  and  for  the  low-frequency 
dielectric  constant  to  be  proportional  to  the  square  root 
of  the  H"^  concentration.  This  is  shown  in  Figure  24,  in 
which  are  plotted  the  50-Hz  values  of  log  conductivity 
(circles)  and  log  dielectric  constant  (crosses)  as  func¬ 
tions  of  pH  (the  negative  log  H'*'  ion  concentration)  of  the 
meltwater  of  the  various  snows.  Other  ionic  impurities 
seem  to  have  little  effect.  This  suggests  that  the  low- 
frequency  conductivity  is  the  result  of  proton  hopping. 
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'8.  Dielectric  spectra  for  sample  16  for  indicated  water  contents  (or  temperatures). 


Dielectric  spectra  for  sample  10  for  indicated  water  contents  (or  temperatures). 
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Figure  20.  Dielectric  spectra  for  sample  16  for  indicated  water  contents  (or  temperatures). 
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Figure  24.  Log  iq  low-frequency  conductivity  (left-hand 
scale)andlogiolow-frequency  dielectric  constant(right- 
harui  scale)  for  various  snows  vs  pH  of  the  meltwater. 


Figure  26.  Logjo  percent  water  vs  logio  dielec¬ 
tric  constant  at  5  kHz  for  the  snows  indicated  by 
Roman  numerals. 


Meltwater  Conductivity  (nS/M) 

Figure  25.  Concentrations  of  sodium  arui  chlor¬ 
ide  ions  vs  conductivity  (I  kHz)  of  the  meltwater 
for  various  snows. 

These  data  have  been  corrected  for  temperature.  To 
correct  for  density  is  more  tenuous  because,  as  will  be 
seen,  our  density  effect  measurements  gave  rather  ragged 
curves  and  were  all  made  at  -9‘’C.  When  we  do  attempt 
to  correct  for  density,  the  points  are  moved  as  indicated 
by  the  arrows.  Except  for  snow  V,  the  general  tenden¬ 
cies  remain.  But  the  effect  on  snow  V  is  catastrophic. 
Possibly  grain  size  plays  a  role  here  because  samples  IV 
and  V,  for  which  the  correction  is  largest,  are  the 
samples  with  the  finest  grains.  Exploring  the  questions 
raised  by  Figure  24  may  be  a  fruitful  area  for  future 
research. 

The  relationship  between  low-frequency  conduc¬ 
tance  of  the  snow  and  the  hydrogen  ion  concentration 
helps  to  explain  the  great  difference  between  snow  lA 


(Fig.  14-16)  and  snow  IB  (Fig.  17, 18).  The  very  low, 
low-frequency  conduction  of  the  latter  exposes  struc¬ 
ture  in  the  log  a  vs  logjo  frequency  curve,  which  is 
obscured  by  the  high  low-frequency  conduction  of 
sample  lA.  The  conductivity  of  the  melt  is  dominated  by 
the  ions  of  NaQ  except  at  its  lowest  concentrations 
where  the  ion  may  play  a  role.  This  can  be  seen  in 
Table  1  and  Figure  25. 

The  water-content  experiments  all  produced  ordered 
sets  of  data.  Even  experiments  on  different  samples  of 
the  same  snow  gave  results  that  fit  together  rather  well. 
In  all,  both  the  conductivity  and  dielectric  constant 
increased  as  water  content  increased.  Thus  it  is  tempting 
to  take  vertical  cuts  of  these  data  to  see  if,  at  a  single  fre¬ 
quency,  a  relation  can  be  found  between  either  param¬ 
eter  and  the  water  content.  Figure  26  is  such  a  cut.  It 
shows  the  variation  with  water  content  of  the  dielectric 
constant  at  5  kHz  on  a  log-log  plot.  For  a  given  snow,  the 
data  fall  in  a  well-order‘»i  sequence  on  a  straight  line. 
Where  we  report  data  on  more  than  one  sample  of  the 
same  snow,  different  symbols  are  used  for  each  sample 
but  the  points  fall  on  nearly  the  same  lines.  This  consis¬ 
tency  builds  confidence  in  our  techniques.  However,  no 
two  different  snows  lie  on  the  same  line.  We  have  found 
no  slice  of  this  kind  that  appears  to  be  useful  for  deter¬ 
mining  water  content.  Although  there  is  an  easily  mea¬ 
surable  variation  of  both  the  dielectric  constant  and 
conductivity  spectra  with  water  content,  no  relation  has 
been  found  that  works  well  with  all  snows.  There  appear 
to  be  too  many  variables  to  sort  out. 

We  have  tried  several  strategies  to  discover  a  way  to 
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Table  3.  Polynomial  fit  to  data  for  samples  10-16. 


Experiment 

number 

Snow 

number 

Meltwater 

(a) 

Melt 

(pH) 

Density 

(StbnL) 

Temp  (T) 
or  %  H-yO 

Polynomial  fit 

af-o^+bfo+cf<^^ 

ABC 

Average 
percent 
relative  error 

10 

V 

68 

4.8 

0.182 

7  =  ^.5 

674.07 

-7.585 

0.0315 

13.8 

10 

V 

68 

4.8 

0.182 

7  = -2 

532.39 

-3.651 

0.0170 

1.8 

10 

V 

68 

4.8 

0.8 

675.41 

-5.683 

0.0234 

4.2 

10 

V 

68 

4.8 

1.7 

609.69 

-3.725 

0.0161 

1.0 

10 

V 

68 

4.8 

2.7 

640.11 

-3.756 

0.0157 

0.7 

10 

V 

68 

4.8 

3.7 

718.70 

^.366 

0.0167 

2.2 

11 

IV 

42 

5 

0.265 

7  =  -5 

370.94 

1.292 

-0.0005 

5.2 

11 

IV 

42 

5 

0.265 

7  =  -5 

234.33 

4.388 

-0.0115 

0.2 

11 

IV 

42 

5 

0.8 

443.06 

3.377 

-0.0111 

3.9 

11 

IV 

42 

5 

2.5 

512.71 

3.279 

-0.01 14 

5.6 

11 

IV 

42 

5 

2.9 

511.56 

3.849 

-0.0151 

6.1 

11 

IV 

42 

5 

3.6 

542.40 

3.970 

-0.0149 

6.1 

12 

II 

220 

4.27 

0.36 

7  =  -2 

634.04 

4.019 

-0.0055 

2.8 

12 

II 

220 

4.27 

0.36 

7  =  -^ 

508.96 

8.705 

-0.0241 

1.8 

12 

II 

220 

4.27 

0.7 

1263.95 

7.472 

-0.0327 

4.6 

12 

II 

220 

4.27 

2.1 

1743.51 

6.002 

-0.0348 

6.0 

12 

II 

220 

4.27 

3.3 

2269.86 

2.865 

-0.0302 

6.3 

12 

II 

220 

4.27 

4.0 

2331.67 

3.313 

-0.0359 

6.1 

13 

Ib 

294 

6.05 

0.423 

7= -3.5 

83.21 

10.048 

-0.0239 

1.3 

13 

Ib 

294 

6.05 

0.423 

II 

) 

93.32 

11.561 

-0.0298 

0.6 

13 

Ib 

294 

6.05 

0.5 

48.83 

16.842 

-0.0505 

2.1 

13 

Ib 

294 

6.05 

1.3 

39.39 

21.434 

-0.0691 

1.5 

13 

Ib 

294 

6.05 

1.8 

45.76 

25.563 

-0.0856 

1.3 

13 

Ib 

294 

6.05 

2.3 

40.62 

28.591 

-0.0985 

1.2 

14 

III 

24 

5.36 

0.391 

7  =  -5 

240.13 

6.089 

-0.0178 

0.8 

14 

III 

24 

5.36 

0.391 

7  = -5 

207.33 

8.457 

-0.0263 

1.9 

14 

III 

24 

5.36 

0.9 

355.49 

10.704 

-0.0397 

6.6 

14 

III 

24 

5.36 

1.6 

423.59 

11.452 

-0.0440 

6.6 

14 

III 

24 

5.36 

1.9 

457.73 

11.788 

-0.0461 

6.7 

14 

in 

24 

5.36 

2.1 

464.52 

12.482 

-0.0493 

6.5 

15 

II 

220 

4.27 

0.372 

7  =  -7 

1134.67 

-9.124 

0.0424 

9.5 

15 

II 

220 

4.27 

0.372 

7= -5 

613.35 

6.700 

-0.0188 

3.8 

15 

II 

220 

4.27 

0.3 

1141.45 

4.213 

-0.0172 

4.2 

15 

II 

220 

4.27 

1.2 

1563.68 

-0.001 

-0.0059 

3.2 

15 

II 

220 

4.27 

2.3 

2309.29 

-7.400 

0.0134 

2.7 

15 

11 

220 

4.27 

3.2 

2933.29 

-14.816 

0.0344 

0.4 

16 

Ib 

294 

6.05 

0.422 

7  =  -8.5 

19.59 

10.429 

-0.0239 

3.1 

16 

Ib 

294 

6.05 

0.422 

11 

86.51 

12.178 

-0.0309 

0.9 

16 

Ib 

294 

6.05 

0.5 

32.84 

18.490 

-0.0556 

2.1 

16 

Ib 

294 

6.05 

0.7 

31.82 

19.470 

-0.0592 

2.2 

16 

Ib 

294 

6.05 

0.9 

29.27 

20.694 

-0.0641 

2.2 

16 

Ib 

294 

6.05 

I.I 

25.21 

22.366 

-0.0708 

2.6 

use  our  data  to  determine  water  content  of  snow.  De¬ 
spite  clear  trends  of  the  dielectric  data  with  water  con¬ 
tent,  they  are  not  free  enough  from  variations  due  to 
other  parameters  to  allow  us  to  work  backward  from 
dielectric  measurements  to  water  content.  One  scheme 
was  to  fit  our  dielectric  constant  data  with  a  polynomial 
of  as  low  an  order  as  we  could  in  the  hope  that  the 
coefficients  would  serve  to  catalog  our  samples.  If  there 
were  orderly  variations  in  these  that  could  be  related  to 
water  content,  then  we  might  have  a  useful  tool.  We 
found  that  we  could  fit  all  of  our  data  for  samples  10-16 
using  an  equation,  Af^'^  +  fi  +  C/°'^,  where  /  is  the 


frequency  in  hertz  and  A,  B  and  C  are  the  parameters 
(see  Table  3).  We  could  fit  the  data  for  samples  1 , 2  and 
3  with  an  equation  Af~'+  +  C.  However,  we 

were  unable  to  find  a  three-parameter  polynomial  form 
which  would  fit  all  of  our  data.  Note  that  the  snow  used 
in  samples  1 , 2  and  3  was  of  type  lA,  having  very  large 
grains  and  large  conductivity. 

Another  approach  was  to  use  the  data  describing  the 
form  of  the  various  Cole  plots,  that  is,  the  center  fre¬ 
quencies  of  the  major  groups  of  dispersions,  the  strengths 
of  these  groups  as  indicated  by  the  diameters  of  the  Cole 
circles,  and  the  distance  below  K"  =  0of  the  centers  of 
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Table  4.  Cole  plot  parameters  for  the  various  samples. 


Low  dispersion  Hinh  dispersion 


Sample 

no. 

Snow 

type 

Condition 

(ohm-m) 

pH 

Temp. 

or%W 

/ 

Hz 

S 

pF 

H 

pF 

/ 

Hz 

s 

pF 

H 

pF 

1 

la 

1.69E-03 

4.4 

T  =  -6 

60 

160 

40 

150 

220 

15 

1 

la 

1.69E-4)3 

4.4 

3.5 

100 

1200 

520 

500 

420 

no 

1 

la 

1.69E-03 

4.4 

24.0 

70 

3600 

1000 

200 

3600 

700 

1 

la 

1.69E-03 

4.4 

38.0 

110 

4600 

350 

2 

la 

1.69E-03 

4.4 

T=  -4.5 

100 

230 

25 

2 

la 

1.69E-03 

4.4 

9.7 

120 

1600 
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Notes: 

S  =  Strength  of  the  dispersion  (diameter  of  Cole  circle) 
H  =  Displacement  of  the  center  of  the  circle  below  1C  = 
f  =  Frequency  corresponding  to  center  of  circle. 

these  circles  (which  is  an  indication  of  the  breadth  of  the 
distribution  of  the  dispersions).  In  the  case  of  two 
groups  of  dispersions,  this  would  require  a  total  of  six 
parameters  per  sample.  These  data  are  given  in  Table  4. 

No  strong  dependence  on  water  content  is  evident,  al¬ 
though  there  does  seem  to  be  a  clear  tendency  for  the 
strength  of  the  dispersions  to  increase  with  water  con¬ 
tent  for  those  samples  having  no  clear  high-frequency 
dispersion.  Again  we  find  a  profound  difference  be¬ 
tween  snow  lA  and  all  of  the  others. 


THE  EFFECT  OF  DENSITY 

A  series  of  experiments  in  which  samples  of  dry 
snow  at  -9°C  were  subjected  to  unconfined  compres¬ 
sion  was  undertaken  to  determine  how  sample  density 
affected  the  dielectric  properties.  A  guarded  ring  cell, 
one  electrode  of  which  was  driven  by  a  micrometer,  was 
constructed  of  stainless  steel  (Fig.  27).  The  cell  was 
placed  on  its  side  in  an  aluminum  box  which  was  loosely 
filled  with  snow.  Bridge  measurements  were  then  made 


27 


Snow  loaded  here 


Bottom  Plate 

|< - 1.5  in. - 


Figure  27.  Micrometer  cell  used  for  variable  density  experiments. 
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at  50  Hz,  500  Hz,  5  kHz  and  50  kHz.  The  micrometer  results  are  shown  in  Figures  28-36.  Figures  28-32  show 

was  turned  to  somewhat  compress  the  sample  and  a  new  frequency  spectra  for  conductivity  and  dielectric  con- 

frequency  spectrum  was  run.  This  process  was  repeated  stant  with  density  as  a  parameter.  Both  increase  some- 

several  times  to  produce  data  over  a  range  of  densities.  what  regularly  with  density  in  a  manner  roughly  propor- 

In  all,  samples  of  four  different  snows  were  run  in  this  tional  to  the  initial  magnitude.  The  scatter  in  the  data 

manner  with  two  different  samples  of  one  snow.  The  from  sample  to  sample  is  emphasized  when  the  conduc- 


15  2  25  3  35  4  455  15  2  25  3  35  4  45  5 


UKflCQUENCY  LOG  reOJENCY 

Figure  28.  Conductivity,  a,  and  dielectric  constant,  K,  as  functions  of 
frequency  for  sample  5.  Density  in  g/cn^  is  shown  as  a  parameter. 
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Figure  29.  Conductivity  and  dielectric  constant  as functions  of  frequency 
for  sample  6.  Density  in  g/cm^  is  shown  as  a  parameter. 
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Figure  30.  Conductivity  and  dielectric  constant  as  functions  of  frequency 
for  sample  7.  Density  in  g/cm^  is  shown  as  a  parameter. 
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Figure  31 .  Conductivity  and  dielectric  constant  as functions  of  frequency 
for  sample  8.  Density  in  g/cn^  is  shown  as  a  parameter. 
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Figure  32 .  Conductivity  and  dielectric  constant  as functions  of  frequency 
for  sample  9.  Density  in  g/cn^  is  shown  as  a  parameter. 


tivity  and  dielectric  constant  are  plotted  against  density 
for  three  given  frequencies  with  the  sample  as  a  para¬ 
meter  (Figures  33-36).  Sample  8  (snow  IB)  is  clearly 
anomalous  as  is  sample  9  (snow  11)  to  a  lesser  extent. 
Each  snow  behaves  in  a  fairly  orderly  way  but  each  is 
somewhat  unique.  To  illustrate  this  more  fully.  Table  S 
gives  the  ratio  of  K  for  maximum  density  to  K  for  mini¬ 
mum  density  and  the  same  for  conductivity.  This  shows 


the  trends  and  the  variability.  At  the  bottom  of  Table  5 
is  given  the  average  for  the  first  four  samples,  which 
have  approximately  the  same  density  range.  This  gives 
an  overall  idea  of  how  density  differences  among  sam¬ 
ples  may  affect  their  dielectric  response. 

Figure  37  is  a  plot  of  the  difference  between  the  con¬ 
ductivity  at  SO  kHz  and  that  at  SO  Hz  vs  density  for  dry 
snow.  This  difference  in  conductivity  is  roughly  pro- 


Figure  33.  Conductivity  and  dielectric  constant  at  50  Hz  vs  density  for  the  samples  indicated 
by  the  number. 
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Figure  34.  Conductivity,  a  (nS/m)  and  dielectric  constant,  K,  at  500  Hz  vs  density  for  the 
samples  indicated  by  the  number. 


Figure  35.  Conductivity  and  dielectric  constant  at  5  Hz  vs  density  for  the  samples  indicated 
by  the  number. 


Figure  36 .  Conductivity  and  dielectric  constant  at  50  kHz  vs  density  for  the  samples  indicated 
by  the  number. 


Figure  37.  Conductivity  difference,  o(\iSI 
m)at50  kHz  minus  that  at  50  Hz  vs  density 
(glcm^)  for  dry  snow.  The  numbers  refer  to 
density  experiments  at  -9  °C  .  The  *  are 
from  samples  10-16  in  the  temperature 
range  from  -2°C  to -5  °C  . 
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portional  to  density,  not  surprising  when  one  considers 
that  it  represents  the  parallel  sum  of  the  conductances 
which  appear  in  all  of  the  Debye  dispersions  making  up 
the  sample.  Because  each  sample  has  a  different  dis¬ 
tribution  of  dispersions,  we  would  expect  each  to  be¬ 
have  differently  in  detail.  Note  that  the  data  converge 
toward  the  origin  as  they  should.  We  were  surprised  that 
a  similar  plot  of  the  dielectric  constant  difference  vs 
density  (Fig.  38)  showed  a  very  broad  scatter  and  only 
a  very  general  tendency  to  increase  with  density. 

However  for  dry  snow  close  to  the  melting  point,  we 
see  in  Figure  39  that  the  logjo  of  dielectric  constant  at  5 
kHz  increases  nearly  proportionally  to  the  logio  of 
density.  Thus  the  dielectric  constant  at  5  kHz  is  roughly 
proportional  to  the  density.  Unfortunately,  the  band  of 
values  in  both  Figures  37  and  39  is  too  broad  to  make 


Figure  38.  Dielectric  constant  at  50  kHz 
minus  that  at  50  Hz  vs  density  for  dry  snow. 
Points  •  are  for  samples  10-16.  The  num¬ 
bers  indicate  samples  used  in  the  density 
experiments. 


Figure  39.  Log  jo  K  (at  5  kHz)  vs  logjo  den¬ 
sity  for  dry  snow  samples  numbered  as 
shown  (-9°C). 


Figure  40.  Log  ;oK  (at  5  kHz)  vs  log jo  den¬ 
sity  for  samples  10-16  having  water  con¬ 
tents  between  1  and  3%.  Where  more  than 
one  X  appears  for  a  given  density,  the  high¬ 
er  the  mark,  the  higher  the  water  content. 

possible  more  than  a  rough  estimate  of  density  from 
measurements  of  o  or  K.  Similar  data  for  samples  hav¬ 
ing  a  small  water  content,  from  1  %  to  3%,  are  shown  in 
Figure  40.  The  same  tendency  is  present;  but  the  values 
of  /IT  are  shifted  upward  by  roughly  50%  from  those  for 
dry  snow. 


CONCLUSIONS  AND  RECOMMENDATIONS 

Detailed  discussion  will  be  found  in  each  appropri¬ 
ate  section.  Here  we  summarize  our  major  observations 
and  comment  on  possibly  useful  further  studies. 

Technique 

The  electrical  heating  procedure  developed  for  in¬ 
crementally  melting  the  snow  in  order  to  vary  the  water 
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content  works  very  well.  By  improving  the  powermeas- 
urement.  particularly  by  improving  the  measurement  of 
phase,  control  of  water  content  can  be  made  more  sensi¬ 
tive.  While  the  accuracy  of  the  present  system  was  ade¬ 
quate  to  our  needs,  greater  precision  might  make  possi¬ 
ble  very  useful  experiments.  For  example,  at  tempera¬ 
tures  close  to  the  melting  point,  the  surface  properties  of 
ice  are  believed  to  change  markedly.  The  mobility  of 
protons  near  the  surface  is  greatly  enhanced  and  the 
material  is  often  said  to  become  “liquid  like.”  One 
would  like  to  study  this  phenomenon  in  a  system  like 
snow  in  which  the  surface-to-volume  ratio  is  high. 
Well-controlled  electrical  heating  of  the  kind  we  have 
used  should  make  it  possible  to  change  the  temperature 
of  the  whole  sample  uniformly  at  temperatures  close  to 
0°C.  Dielectric  measurements  should  then  be  respon¬ 
sive  to  changes  in  the  surface  properties. 

Our  technique  of  monitoring  the  electrode  separa¬ 
tion  also  makes  it  possible  to  study  the  dimensional 
changes  of  a  system  under  constant  load  as  the  water 
content  changes.  We  suggest  that  interesting  mechani¬ 
cal  information  about  snow,  and  possibly  other  systems, 
may  be  obtained  by  exploring  dimensional  changes  in 
this  way,  particularly  at  temperatures  close  to  the  melt¬ 
ing  point.  We  believe  that  this  technique  could  be  made 
more  sensitive  with  more  sophisticated  apparatus  and 
that  an  effort  at  such  improvement  would  be  worth¬ 
while. 

The  consistency  of  the  data  gives  some  confidence  in 
the  reliability  of  the  technique.  For  example,  the  data  for 
samples  1 , 2  and  3,  all  from  snow  lA,  form  a  consistent 
set  for  1 0  different  water  contents.  Their  spectra  at  zero 
water  content  and  approximately  the  same  temperature 
are  very  nearly  the  same  (see  Fig.  14-16). 

Meltwater 

Room  temperature  analysis  of  the  meltwater  of  the 
various  snows  shows  a  wide  range  of  impurities,  C*  and 
SO4  being  the  most  abundant  ions.  The  details  are  given 
in  Table  1.  Figure  25  shows  that  the  conductivity  is 
dominated  by  NaCl  except  at  the  lowest  values.  It  may 
also  be  of  interest  that  the  total  anions  minus  the  total 
cations  is  the  most  positive  for  the  lowest  pH  and  most 
negative  for  the  two  highest  values  of  pH.  We  believe 
that  a  continuing  program  of  routine  collection  and 
analysis  of  fresh  snow  might  be  useful.  Much  analysis 
of  this  kind  is  done  in  connection  with  acid  rain  studies, 
but  in  the  work  with  which  we  are  familiar,  rain  and 
snow  are  sometimes  collected  together  and  periodically 
analyzed. 

Dry  snow 

For  dry  snow  (at  about  -9°C),  the  low-frequency 
conductivity  is  proportional  to  the  hydrogen-ion  con¬ 


centration  and  seems  quite  insensitive  to  other  ions  (see 
Fig.  24).  Thus  proton  hopping  seems  to  be  the  dominant 
process.  If  this  were  a  surface  conductivity,  we  would 
expect  a  large  effect  of  grain  size,  of  which  we  find  a  hint 
in  the  discussion  of  Figure  24.  At  -9°C,  we  would  not 
expect  large  surface  conductivity  for  ice,  but  it  might  be 
important  as  the  melting  point  was  approached.  As 
noted  above,  our  mechanism  for  heating  would  seem 
well  suited  to  such  studies. 

It  was  shown  in  Figure  37  that  the  difference  between 
the  high-frequency  conductivity  and  the  low-frequency 
conductivity  was  roughly  proportional  to  density  but 
that  the  proportionality  was  different  for  different  sam¬ 
ples.  No  such  relation  seemed  to  hold  for  the  differences 
in  dielectric  constant.  But  there  was  a  proportionality 
for  the  dielectric  constant  measured  at  5  kHz. 

Pressure  of  itself  seems  to  have  little  or  no  effect  on 
the  dielectric  properties  of  dry  snow.  This  would  indi¬ 
cate  that  the  contact  pressure  between  grains  plays  only 
a  small  role  in  determining  the  electrical  properties  of 
dry  snow.  (This  tentative  result  should  be  followed  up 
by  more  extensive  experiments.)  However,  the  change 
in  density  resulting  from  pressure  is  very  important. 

Although  based  on  only  five  different  snows,  our 
measurements  suggest  that  the  low-frequency  conduc¬ 
tivity  of  dry  snow  is  proportional  to  the  hydrogen  ion 
concentration  and  that  the  low-frequency  dielectric 
constant  is  proportional  to  the  square  root  of  the  hydro¬ 
gen  ion  concentration.  To  explore  this  further  would 
seem  to  be  a  fairly  simple  and  useful  project.  If  proton 
hopping  is  the  dominant  process  for  low-frequency 
conduction  in  dry  snow,  activation  energy  measure¬ 
ments  in  this  region  might  provide  some  confirmation 
and  useful  information  relating  to  bulK  and  surface 
processes.  This  too  should  be  a  straightforward  experi¬ 
ment. 

Wet  snow 

Mechanical 

Because  of  the  sensitivity  of  our  apparatus  for  mea¬ 
suring  the  thickness  of  the  sample,  we  were  able  to 
monitor  the  dimensional  change  as  the  water  content 
increased.  There  was  a  steep  change  at  the  start  and  then 
a  rather  uniform  change  as  the  water  content  increased. 
The  steep  change  at  the  start  we  attribute  to  capillarity 
and  possibly  to  liquid-like  surface  states  as  the  tempera¬ 
ture  approaches  very  close  to  zero.  Our  technique  seems 
well  suited  to  studying  this  process  in  more  detail. 

Electrical 

We  have  developed  a  body  of  data  for  five  different 
snows  which  is  as  detailed  and  complete  as  time  and 
facilities  permitted.  It  is  presented  in  a  form  making  it 
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availaDle  for  further  analysis.  The  dominant  properties 
affecting  the  dielectric  spectra  are  water  content,  densi¬ 
ty  and  pH  of  the  melt  water.  Grain  size  seems  to  be  im¬ 
portant,  but  we  do  not  have  enough  detailed  information 
to  be  sure.  Ionic  impurities  do  not  appear  to  play  a  strong 
role  in  dry  snow  but  they  seem  to  make  a  major  con¬ 
tribution  to  the  meltwater  conduction  when  their  abun¬ 
dance  is  high. 

Experiments  on  dry  snow  indicated  that  the  inter¬ 
granular  pressure  was  not  important  in  determining  the 
dielectric  properties.  It  would  be  useful  to  explore  this 
for  wet  snow. 

The  dielectric  properties  vary  strongly  with  water 
content.  Both  conductivity  and  dielectric  constant  in¬ 
crease  but  not  in  the  same  manner  from  sample  to  sam¬ 
ple.  This  is  not  surprising  considering  the  complexity  of 
the  system.  We  have  explored  two  different  approaches 
to  searching  for  useful  regularities  in  it.  In  one,  we  at¬ 
tempted  to  fit  all  of  our  data  with  a  three-parameter 
power  law  polynomial,  each  sartiple  to  be  characterized 
by  different  coefficients.  (Our  thought  was  that  if  there 
were  three  important  properties,  they  might  be  sepa¬ 
rable  in  this  way.)  We  were  able  to  fit  much,  but  not  all, 
of  our  data  with  one  equation.  In  the  other  approach  we 
analyzed  our  data  in  terms  of  Cole  plot  parameters 
(Table  4).  In  neither  did  we  establish  a  method  by  which 
these  data  might  be  used  to  evaluate  the  water  content  of 
an  unknown  sample.  However,  we  think  it  might  be 
interesting  to  try  to  deconvolute  these  spectra  by  assum¬ 
ing  that  water  content,  density,  grain  size  and  possibly 
hydrogen  ion  concentration  each  affect  the  spectra  in  its 
own  characteristic  way  and  that  they  combine  in  a  linear 
way  to  produce  the  measured  spectra. 
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APPENDIX  A.  DATA  TABLE 


Sample 

no. 

Density 

(glcm^) 

Temp  (’C) 
or  %  H2O 

Logy. 

50  Hz 

LogK 
500  Hz 

LogK 

5  kHz 

LogK 

50  kHz 

Logo 

50  Hz 

Logo 
500  Hz 

Logo 

5  kHz 

Logo 
50  kHz 

I 

0.42 

r=-6 

2.473 

1.620 

1.122 

0.802 

-5.539 

-5.344 

-5.161 

-4.824 

1 A 

I 

0.42 

3.5% 

2.835 

2.223 

1.476 

0.455 

-♦.695 

-4.636 

-4.481 

-4.343 

1 

24.0% 

3.698 

2.m 

1.866 

1.062 

-3.965 

-3.884 

-3.784 

-3.736 

1 

38.0% 

3.821 

2.932 

1.893 

1.238 

-3.837 

-3.711 

-3.617 

-3.568 

2 

0.42 

r  =  -4.5 

2.380 

1.538 

1.154 

0.807 

-5.547 

-5.377 

-5.191 

-4.801 

lA 

2 

9.7% 

3.332 

2.623 

1.796 

0.785 

-♦.348 

-4.272 

-4.127 

-4.002 

2 

28.6% 

3.811 

2.987 

1.982 

I.I18 

-3.814 

-3.725 

-3.637 

-3.598 

3 

0.42 

r  =  -4.5 

2.378 

1.678 

1.171 

0.821 

-5.588 

-5.393 

-5.158 

-4.806 

1 A 

3 

r=-o.2 

2.240 

1.641 

1.234 

0.831 

-5.414 

-5.317 

-5.123 

-4.753 

3 

2.5% 

2.881 

2.203 

1.503 

0.777 

-4.690 

-4.614 

-4.472 

-4.310 

3 

5.7% 

3.068 

2.409 

1.634 

0.840 

-4.531 

-4.454 

-4.299 

-4.172 

4 

0.264 

r=-9 

2.029 

1.364 

0.954 

0.689 

-5.585 

-5.520 

-5.397 

-5.032 

4 

0.317 

r=-9 

2.460 

1.742 

I.L0O 

0.737 

-5.078 

-4.990 

-4.909 

-4.699 

II 

4 

0.366 

T  =  -9 

2.635 

1.996 

1.279 

0.823 

-4.806 

-4.730 

-4.645 

-4.465 

4 

0.396 

T  =  -9 

2.808 

2.063 

1.325 

0.867 

^.739 

-4.643 

-4.568 

-4.395 

4 

0.432 

T  =  -9 

2.887 

2.134 

1.385 

0.900 

-4.637 

-4.548 

-4.414 

-4.310 

5 

0.228 

r=-9 

1.457 

0.885 

0.595 

0.320 

-7.106 

-6.685 

-6.222 

-6.632 

5 

0.337 

r=-9 

1.927 

1.282 

0.798 

0.427 

-6.260 

-5.987 

-5.692 

-5.286 

III 

5 

0.405 

T  =  -9 

2.222 

1.517 

0.911 

0.465 

-5.917 

-5.666 

-5.407 

-5.124 

5 

0.476 

T  =  -9 

2.319 

1.626 

0.996 

0.442 

-5.770 

-5.535 

-5.282 

-5.037 

5 

0.543 

T  =  -9 

2.356 

1.692 

1.077 

0.469 

-5.681 

-5.456 

-5.189 

^.951 

6 

0.225 

T  =  -9 

1.383 

0.929 

0.667 

0.376 

-6.972 

-6.664 

-6.169 

-5.483 

6 

0.246 

T  =  -9 

1.330 

0.990 

0.710 

0.394 

-6.630 

-6.452 

-6.032 

-5.443 

III 

6 

0.281 

T  =  -9 

1.450 

1.099 

0.790 

0.407 

-6.380 

-6.230 

-5.864 

-5.314 

6 

0.328 

T  =  -9 

1.607 

1.216 

0.870 

0.476 

-6.156 

-6.027 

-5.696 

-5.215 

6 

0.463 

T  =  -9 

2.034 

1.565 

1.075 

0.511 

-5.621 

-5.510 

-5.241 

-4.971 

7 

0.111 

T^-9 

1.343 

0.901 

0.524 

0.228 

-7.IV5 

-6.632 

-6.231 

-5.657 

7 

0.166 

T  =  -9 

1.325 

0.859 

0.585 

0.317 

-6.004 

-5.959 

-5.833 

-5.471 

IV 

7 

0.248 

T  =  -9 

1.580 

1.001 

0.651 

0.404 

-5.611 

-5.573 

-5.501 

-5.249 

7 

0.497 

T  =  -9 

2.213 

1.628 

0.915 

0.485 

-5.013 

-4.965 

-4.881 

-4.792 

7 

0.662 

T  =  -9 

2.137 

1.649 

0.987 

0.462 

-4.999 

-4.952 

-4.846 

-4.751 

8 

0.309 

T  =  -9 

2.068 

1.344 

0.951 

0.665 

-5.581 

-5.485 

-5.344 

-5.014 

8 

0.355 

T  =  -9 

2.029 

1.280 

0.903 

0.687 

-5.068 

-5.039 

-5.001 

-4.820 

li 

8 

0.440 

T  =  -9 

2.018 

1.403 

0.992 

0.678 

-4.656 

-4.647 

-4.630 

-4.527 

8 

0.500 

T  =  -9 

2.249 

1.572 

l.(»l 

0.680 

-4.513 

-4.498 

-4.478 

-4.384 

8 

0.564 

T  =  -9 

2.219 

1.551 

1.101 

0.673 

-4.452 

-4.437 

-4.413 

-4.326 

9 

0.423 

r  =  -9 

1.087 

0.967 

0.910 

0.711 

-7.410 

-7.128 

-6.436 

-5.159 

9 

0.441 

T  =  -9 

1.230 

1.098 

1.026 

0.733 

-7.091 

-6.884 

-6.181 

-5.017 

IB 

9 

0.470 

T=-9 

1.351 

1.209 

1. 122 

0.758 

-6.938 

-6.734 

-6.006 

-4.912 

9 

0.504 

T  =  -9 

1.477 

1.326 

1.232 

0.791 

-6.799 

-6.589 

-5.850 

-4.799 

9 

0.564 

T  =  -9 

1.539 

1.379 

1.276 

0.798 

-6.708 

-6.494 

-5.778 

^.740 

10 

0.182 

r  =  -4.5 

2.001 

1.250 

0.674 

0.391 

-5.852 

-5.704 

-5.564 

-5.289 

10 

0.182 

r^-0.2 

1.865 

1.297 

0.712 

0.402 

-5.697 

-5.590 

-5.454 

-5.221 

y 

10 

0.8% 

1.975 

1.363 

0.757 

0.408 

-5.789 

-5.649 

-5.483 

-5.226 

10 

1.7% 

1.922 

1.370 

0.784 

0.415 

-5.794 

-5.656 

-5.474 

-5.208 

10 

2.7% 

1.943 

1.396 

0.809 

0.417 

-5.788 

-5.649 

-5.453 

-5.191 

10 

3.7% 

2.000 

1.431 

0.851 

0.410 

-5.758 

-5.615 

-5.410 

-5.166 

11 

0.265 

r=-5 

1.761 

1.206 

0.825 

0.452 

-6.074 

-5.932 

-5.666 

-5.215 

11 

0.265 

r=-5 

1.572 

1.166 

0.838 

0.457 

-5.731 

-5.682 

-5.522 

-5.151 

IV 

11 

0-8% 

1.793 

1.392 

0.939 

0.459 

-5.915 

-5.789 

-5.483 

-5.113 

11 

2.5% 

1.841 

1.460 

0.975 

0.480 

-5.901 

-5.769 

-5.443 

-5.084 

11 

2.9% 

1.839 

1.471 

0.987 

0.441 

-5.896 

-5.762 

-5.425 

-5.084 

11 

3.6% 

1.864 

1.495 

1.011 

0.488 

-5.872 

-5.734 

-5.403 

-5.075 

37 


Sample 

no. 

Density 

(g/cm^) 

Temp(°C) 
or  %  H2O 

LogK 
50  Hz 

UgK 
500  Hz 

LogK 

5  kHz 

LogK 

50  kHz 

Logo 

50  Hz 

Logo 
500  Hz 

Logo 

5  kHz 

Logo 
50  kHz 

i: 

0.36 

T  =  -2 

1.990 

1.484 

1.106 

0.749 

-5.273 

-5.239 

-5.120 

-1.814 

12 

0.36 

T^4 

1.893 

1.505 

1.149 

0.747 

-5.329 

-5.284 

-5.137 

-4.807 

II 

12 

0.7% 

2.239 

1.839 

1.35 1 

0.765 

-5.209 

-5.120 

-4.907 

-4.618 

12 

2.1% 

2.363 

1.968 

1.435 

0.781 

-5.041 

^.968 

-4.764 

-1.513 

12 

3.3% 

2.470 

2.066 

1.499 

0.799 

^.916 

-4.843 

-4.651 

-4.428 

12 

4.0% 

2.483 

2.078 

1.512 

0.759 

^.854 

-4.7S6 

-1.608 

-4.401 

13 

0.423 

7  =  -3.5 

1.348 

1.112 

0.981 

0.706 

-6.950 

-6.646 

-6.109 

-5.080 

13 

0.423 

7  =  -25 

1.395 

1.174 

1.033 

0.725 

-6.947 

-6.635 

-6.052 

-5.023 

IB 

13 

0.5% 

1.391 

1.239 

1.147 

0.760 

-7.192 

-6.865 

-5.996 

-4.909 

13 

1.3% 

1.439 

1.327 

1.234 

0.790 

-7.198 

-6.848 

-5.876 

-4.822 

13 

1.8% 

1.511 

1.402 

1.305 

0.822 

-7.159 

-6.800 

-5.776 

-4.747 

13 

2.3% 

1.540 

1.443 

1.347 

0.830 

-7.164 

-6.787 

-5.723 

-4.708 

14 

0.391 

7  =  -5 

1.607 

1.209 

0.917 

0.503 

-6.764 

-6.395 

-5.829 

-5.195 

14 

0.391 

7  = -5 

1.560 

1.249 

0.976 

0.545 

-6.327 

-6.174 

-5.734 

-5.128 

14 

0.9% 

1.730 

1.460 

l.IOI 

0.536 

-6.679 

-6.270 

-5.586 

-5.044 

14 

1.6% 

1.800 

1.518 

1.145 

0.546 

-6.639 

-6.211 

-5.517 

-5.002 

14 

1.9% 

1.829 

1.546 

1.165 

0.549 

-6.632 

-6.183 

-5.489 

-4.983 

14 

2.1% 

1.840 

1.557 

I.I8I 

0.549 

-6.621 

-6.176 

-5.471 

-4.975 

15 

0.372 

II 

2.226 

1.546 

1.028 

0.734 

-5.636 

-5.487 

-5.301 

-1.941 

15 

0.372 

7  =  -0.5 

1.943 

1.559 

1.140 

0.720 

-5.020 

-4.993 

-4.900 

-1.696 

15 

0.3% 

2.1gl 

1.774 

1.271 

0.739 

-4.975 

-4.938 

-1.816 

-1.602 

15 

1.2% 

2.325 

1.870 

1.328 

0.755 

-4.951 

^.906 

-4.773 

-4.553 

15 

2.3% 

2.488 

2.006 

I.4I0 

0.774 

^.845 

-1.791 

-4.651 

-4.456 

15 

3.2% 

2.600 

2.072 

1.463 

0.779 

-4.774 

-4.722 

-4.584 

-1.392 

16 

0.422 

7  =  -8.5 

1.149 

1.013 

0.957 

0.714 

-1.399 

-7.107 

-6.351 

-5.102 

16 

0.422 

7  =  -5 

1.375 

1.192 

1.049 

0.752 

-6.638 

-6.447 

-5.988 

-1.998 

16 

0.5% 

1.379 

1.259 

1.178 

0.792 

-7.175 

-6.873 

-5.985 

-1.878 

16 

0.7% 

1.395 

1.278 

1.198 

0.804 

-7.194 

-6.885 

-5.962 

-1.855 

16 

0.9% 

1.411 

1.300 

I.22I 

0.812 

-7.222 

-6.896 

-5.934 

-4.832 

16 

1.1% 

1,433 

1.325 

1.250 

0.823 

-7.222 

-6.901 

-5.895 

-4.802 
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